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Symbols 


A 

nuclear mass number 

b 

impact parameter, fm 

NRL 

Naval Research Laboratory 

Z 

nuclear charge number 

AA 

total number of abraded and ablated nucleons 

^abl 

number of ablated nucleons 

^abr 

number of abraded nucleons 

O 

cross section, mb 

(T 0 , n, 7] 

Subscripts: 

parameters in Silberberg-Tsao theory (eq. (10)) 

F 

fragment 

FSI 

frictional spectator interaction 

nuc 

nuclear 

P 

projectile 

T 

target 
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Abstract 

Cross. S ectwn predictions with semiempirical nuclear fragmentation 
odels from the Langley Research Center and the Naval Research Labo 
ratory are compared with experimental data for the breakup of relativistic 
iron and argon projectile nuclei in various targets. Both these models 
are commonly used to provide fragmentation cross-section inpuTs into 

olemii ‘""I?? ^ ,0T ‘ hieU ' nS and exp ° mr * 

experimental data ' ‘ W “” a «reemenl will, the 


Introduction 


In the approaching era of career astronauts and 
space workers who will man Space Station Freedom , 
establish lunar bases, and explore the solar system, 
concern is mounting over possible deleterious effects 
to crews from the heavy ion component of solar and 
galactic cosmic rays (refs. 1 and 2). To properly as- 
sess these risks, knowledge of cosmic ray interaction 
and transport in bulk matter is required to accu- 
rately determine shielding requirements and to ad- 
equately assess radiobiological damage to the astro- 
nauts. A major source of uncertainty in these risk 
assessments is the input fragmentation cross-section 
data base (ref. 3). At present, the experimental data 
base is inadequate, and accurate theories of nuclear 
fragmentation are hampered by the paucity of experi- 
mental data. Two nuclear fragmentation models cur- 
rently used for galactic cosmic ray shielding studies 
m-e semiempirical formalisms developed at the Naval 
Research Laboratory (refs. 4, 5, and 6) and at the 
Langley Research Center (ref. 7). 

The NRL model involves extrapolations to heavy 
targets of a modified form of a parameterization orig- 
mally developed by Rudstam for hydrogen targets 
fret. 8). Numerous adjustable parameters have been 
chosen by comparisons with available experimental 
data. The Langley model is based upon a two-step 
abrasion-ablation collision formalism. It has one ad- 
justable parameter, a second-order correction to the 
excitation energy used as input into the ablation 
stage of the reaction. 

In the present work, cross-section predictions 
from each semiempirical model are made and com- 
pared with available data from recent experiments 
using iron (ref. 9) and argon beams (ref. 10). Com- 
parisons with earlier measurements (ref. 11) for iron 
beams at energies different from those used in ref- 
erence 9 are also made. The agreement between 
model predictions and experimental measurements is 
assessed by analyzing the distribution of cross-section 


Semiempirical Models 

Formulation of Langley Research Center 

Model 

In the Langley semiempirical model, the classi- 
cal, geometric abrasion-ablation model of Bowman 
Swiatecki, and Tsang (ref. 12) is modified to include 
frictional spectator interactions through the use of 
higher order corrections to the abraded prefragment 
excitation energies. In this method, the nuclear frag- 
mentation cross sections are given by 


<rnuc(Z F , A f ) = F x exp (~R\Z F - SA F + TA 2 F \ 3 / 2> j a (AA) 

( 1 ) 

where according to Rudstam (ref. 8), R = 11.8A" 


- 0.45 


S 0.486, T 3.8 x 10 4 , and F\ is a normalizing 
factor such that 


A f ) = <r(A A) (2) 

z F 


which ensures charge and mass conservation. The 
Rudstam formula for o(AA) is not used because the 
AA dependence is too simple and breaks down for 
eavy targets. Instead, the cross section for removal 
of A A nucleons is estimated by using 

a(AA) = 7T&2 - 7 xb\ (3) 

where b 2 is the impact parameter at which A abr nu- 
cleons are abraded by the collision and A abl nucleons 
are ablated m the subsequent prefragment deexcita- 
tion, such that 


A a br(^2) + A ab ](6 2 ) = A A — ^ 
and similarly for b\ 


^abr(^l) + A ab] (f>i) = A A + i 


( 4 ) 

( 5 ) 


The number of abraded nucleons is estimated from 
the geometric overlap volume and the mean-free path 
in nuclear matter A as 


A abr = FA P 


1 - 0.5 exp 



0.5 exp 



where F is the fraction of the volume in the geomet- 
ric overlap region between the colliding nuclei and 
C P and Cp are the maximum chord lengths of the 
intersecting surfaces in the projectile (P) and target 
(T). Expressions for F given elsewhere (ref. 7) dit- 
fer because of the relative sizes of the colliding nuclei 
and the nature of the collision (central versus pe- 
ripheral). The number of ablated nucleons A a bi ls 
computed from 

* Es ~h ^FSI 

Aabl ~ 10 MeV 

which assumes that a nucleon is ablated (evaporated) 
for every 10 MeV of excitation energy. In equa- 
tion (7), E s represents excitation energy associated 
with the surface energy contribution from abrasion, 
and Ep si represents the contributions resulting from 
frictional spectator interactions. The only arbitrarily 
adjusted parameter in this model is a second-orc er 
correction to the expression for the surface energy 
term. 

Because the dissociation of projectile and target 
nuclei by their interacting Coulomb fields may be 
important for some heavier nuclei at high energies, 
the electromagnetic dissociation contributions cr em 
must be added to the nuclear fragmentation cross 
section <r n uc to yield the total fragmentation cross 

section , , 

Of — C^nuc T <^em ^ ' 

Methods for estimating cr em have been developed and 
parameterized for use with this fragmentation model 
(refs. 7 and 13). 

Formulation of Naval Research 

Laboratory Model 

The fragmentation cross sections for nucleus- 
nucleus collisions with the NRL model are calculated 
from nucleus-nucleon collisions by 

crp(Ap — Ap) = op{Ap — H)5 c £n£L e l e A (®) 

where <J F (A P - H) is the fragmentation cross section 
for nuclear breakup by hydrogen targets. In equa- 
tion (9), S c is a scaling factor obtained by empir- 
ically fitting nuclear skin thicknesses. The factors 


£n and e A> respectively, represent adjustable 

correction factors for neutron- deficient fragments, for 
light mass products, for single-nucleon stripping, and 
for large A A removal. Parameterized expressions for 
these factors and their appropriate limits of applica- 
bility can be found in references 4, 5, and 6. 

From reference 4, the cross sections for fragmen- 
tation by hydrogen targets are given by 

crp(Ap - H) = cr 0 f(Ap)f{E) exp( — PAT) 

x exp (~R\Z - SAp + TA 2 f \ v ) 

( 10 ) 

Equation (10) is applicable to projectile nuclei with 
mass numbers between 9 and 209 and fragments with 
mass numbers Af between 6 and 200, except or 
peripheral interactions where AA(= Ap — Ap) is 
small. Parameterizations of the various factors in 
equation (10) are given elsewhere (refs. 4 and 5). 

Cross-Section Predictions 

With the Langley and NRL semiempirical mod- 
els elemental production cross sections for iron 
beams at 1.88,4 GeV and 1.55 A GeV fragmenting 
in various targets are presented in tables 1 and 2. 
The experimental data are taken from Westfall et al. 
(ref. 11) and Cummings et al. (ref. 9). From ta- 
bles 1 and 2, generally good agreement exists between 
the Langley model predictions and the experimental 
measurements. The NRL model predictions, how- 
ever, typically overestimate the experimental data, 
especially for heavier mass fragments. Detailed anal- 
yses of the distributions of cross-section differences 
are presented in the next section. 

Recently, Tull reported measurements of frag- 
ment production cross sections for 1.65 A GeV argon 
beams fragmenting in carbon and potassium chlo- 
ride (KC1) targets. (See ref. 10.) Figures 1 and 2 
display predictions of elemental productipn cross sec- 
tions obtained with the Langley and the NRL models 
compared with the measured values of Tull. Unlike 
the previous comparisons involving iron beams, the 
agreement between theory and experiment is good 
for both the Langley and the NRL models. Al- 
though not displayed here, comparisons between the- 
ory and experiment were also made for fragment iso- 
tope production cross sections. Detailed analyses of 
the distributions of both elemental and isotopic cross- 
section differences are presented in the next section. 

Distributions of Cross-Section 
Differences 

Quantitative agreement between theory and 
experiment is evaluated by investigating the 
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Element produced, Z 

Figure 1. Elemental production cross sections for 1.65/1 GeV 
argon beams fragmenting in carbon targets. 



Element produced, Z 

Figure 2. Elemental production cross sections for 1.65/1 GeV 
argon beams fragmenting in KC1 targets. 


distribution of cross-section differences. Deciding 
whether theory and experiment agree or disagree is 
actually a subjective interpretation of the results of 
the evaluation process. For example, in some applica- 
tions, differences of up to 50 percent may be consid- 
ered acceptable. For other applications, differences 
greater than 25 percent may not be acceptable. 

In tables 3 and 4, the target-averaged distribu- 
tions of elemental cross-section differences are tab- 
ulated for each incident beam-energy combination. 
The table entries are the percentage of cross-section 
differences within the experimental uncertainties; the 
percentages outside the error bars but within 10, 25, 
50, and 100 percent; and the percentages which differ 
by more than 100 percent. 


From table 3, for 1.88/ GeV iron beams, 62 per- 
cent of the Langley cross-section predictions fall 
within the experimental uncertainties, 77 percent of 
the predictions fall within 25 percent of the experi- 
mental data, and 95 percent fall within 50 percent of 
the data. None of the Langley cross-section predic- 
tions differ by more than 100 percent from the data. 
For the NRL model, 17 percent of the predictions 
fall within the experimental error, 25 percent within 
25 percent of the data, 47 percent within 50 percent 
of the data, and 15 percent differ from the data by 
more than 100 percent. 

The 1.55,4 GeV iron-beam comparisons, also pre- 
sented in table 3, indicate that both models predict 
5 percent of the cross sections falling within the ex- 
perimental errors. For the Langley model, 24 percent 
fall within 10 percent of the data, 59 percent within 
25 percent of the data, and 95 percent within 50 per- 
cent of the data. None of the Langley predictions dif- 
fer from the data by more than 100 percent. For the 
NRL model, 7 percent of the cross-section predictions 
fall within 10 percent of the data, 14 percent agree 
within 25 percent, 45 percent agree within 50 per- 
cent, and 21 percent differ by more than 100 percent. 
Overall, the Langley model appears to give much bet- 
ter agreement with experiment for these iron beams 
fragmenting in various heavy targets. 

In table 4, results for elemental and isotopic cross- 
section differences are presented for 1.65/ GeV ar- 
gon beams fragmenting in carbon and KC1 targets. 
For the Langley model, 8 percent of the elemental 
cross-section predictions fall within the experimental 
uncertainties, 20 percent are within 10 percent of the 
data, 54 percent are within 25 percent, 81 percent are 
within 50 percent, and 92 percent are within 100 per- 
cent of the experimental data. For the NRL model. 

4 percent of the elemental cross-section predictions 
fall within the experimental uncertainties, 16 per- 
cent are within 10 percent of the data, 58 percent 
are within 25 percent, 73 percent are within 50 per- 
cent, and 92 percent are within 100 percent of the 
experimental data. 

Comparing isotopic cross sections, 35 percent of 
the Langley model predictions fall within the error 
bars, 40 percent are within 25 percent of the data, 

53 percent arc within 50 percent, and 89 percent 
are within 100 percent of the experimental values. 

For the NRL model, 34 percent of the isotopic cross 
sections are within the error bars, 41 percent are 
within 25 percent of the data, 59 percent are within 
50 percent, and 81 percent are within 100 percent 
of the experimental data. Overall, these two models 
appear to yield essentially the same agreement with 
experiment for these argon data. 
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Concluding Remarks 

The cross-section predictions of two semi- 
empirical fragmentation models have been com- 
pared with experimental measurements for relativis- 
tic beams of iron and argon colliding with various 
targets. Overall, the Langley Research Center model 
appears to yield better agreement with these data. 
Incorporating the Langley model into cosmic ray 
transport codes should provide improved accuracy 
in predictions of radiation exposures and concomi- 
tant shield requirements for spacecraft crews. For 
elemental production, the Langley model typically 
predicted cross sections which were within 25 per- 
cent of the experimental values for over 80 percent 
of these cross sections. For isotopic production, the 
Langley model had a 53- percent success rate for pre- 
dicting cross sections within 50 percent of the data 
and an 89-percent success rate for predicting cross 
sections within 100 percent of the data. Further com- 
parisons with experiment require additional experi- 
mental data. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
March 23, 1993 
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Table 1. Element Production Cross Sections for 1.88,4 GeV 
Iron Beams Fragmenting in Various Targets 


Element 

produced 

Cross section, mb 

NRL 

Langley 

Experiment 1 

Carbon target 

Mn 

237 

184 

181 ±27 

Cr 

182 

123 

124 ± 13 

V 

115 

101 

100 ±11 

Ti 

157 

87 

87 ± 11 

Sc 

116 

78 

54 dr 9 

Ca 

111 

71 

78 ± 11 

K 

81 

65 

52 i 7 

Ar 

82 

61 

55 dr 9 

Cl 

55 

57 

53 ± 7 

s 

62 

53 

54 ± 10 

p 

40 

50 

59 ± 10 

Si 

39 

47 

57 ± 10 

A1 

31 

44 

83 ±11 

Sulp] 

iur target 

Mn 

402 

217 

250 ± 22 

Cr 

213 

139 

128 ± 16 

V 

135 

115 

86 ± 12 

Ti 

184 

100 

64 ± 10 

Sc 

136 

90 

91 ± 13 

Ca 

130 

82 

97 ± 14 

K 

95 

76 

55 ±21 

Ar 

96 

71 

74 ± 13 

Cl 

64 

66 

66 ± 14 

s 

72 

62 

74 ± 12 

p 

47 

59 

50 ± 8 

Si 

46 

56 

106 ± 14 

A1 

36 

53 

78 ± 13 

Copper target 

Mn 

648 

266 

219 ± 20 

Cr 

250 

158 

149 ± 16 

V 

158 

132 

121 ± 15 

Ti 

216 

117 

101 ± 14 

Sc 

160 

106 

100 ± 15 

Ca 

153 

98 

98 ± 14 

K 

112 

91 

88 ± 14 

Ar 

112 

86 

95 ± 15 

Cl 

75 

82 

86 ± 13 

S 

85 

78 

56 ± 11 

p 

55 

74 

88 ± 15 

Si 

54 

72 

72 ± 11 

A1 

42 

69 

179 ± 27 


a Data from reference 11. 






Table 2. Element Production Cross Sections for 1.55,4 GeV 
Iron Beams Fragmenting in Various Targets 



Cross section, mb 

Element 

produced 

NRL 

Langley 

Experiment 0 

Carbon target 

Mn 

243 

185 

140.73 ± 3.36 

Cr 

196 

124 

105.33 ± 2.69 

V 

121 

101 

79.32 ± 2.31 

Ti 

162 

87 

75.17 ± 2.23 

Sc 

118 

78 

57.29 ± 1.92 

Ca 

111 

71 

63.37 ± 2.01 

K 

80 

65 

43.62 ± 1.64 

Ar 

79 

60 

47.65 ± 1.72 

Cl 

52 

56 

41.45 ± 1.59 

S 

58 

53 

46.47 ± 1.68 

P 

38 

49 

39.45 ± 1.53 

Si 

36 

47 

50.99 ± 1.75 

A1 

27 

44 

41.23 ± 1.55 

Mg 

29 

42 

45.45 ± 1.62 

Na 

24 

40 

35.83 ± 1.42 

Ne 

25 

37 

44.79 ± 1.59 

Aluminum target 

Mn 

359 

208 

174.04 ± 4.46 

Cr 

223 

137 

127.60 ± 3.23 

V 

137 

113 

91.05 ± 2.70 

Ti 

184 

98 

84.12 ±2.58 

Sc 

134 

87 

73.41 ± 2.40 

Ca 

126 

79 

68.92 ± 2.31 

K 

91 

74 

52.89 ±2.01 

Ar 

90 

68 

52.72 ±2.01 

Cl 

59 

64 

45.24 ± 1.85 

S 

66 

60 

52.27 ± 1.98 

P 

43 

57 

43.47 ± 1.80 

Si 

41 

54 

58.21 ± 2.08 

A1 

31 

51 

45.37 ± 1.82 

Mg 

33 

49 

51.76 ± 1.94 

Na 

27 

46 

45.23 ± 1.81 

Ne 

29 

44 

49.11 ± 1.88 


a Data from reference 9. 




Table 2. Concluded 


Cross section, mb 


Element 

produced 

NRL 

Langley 

Experiment 0 

Copper target 

Mn 

670 

263 

238.96 ± 

6.78 

Cr 

270 

159 

147.44 ± 

3.73 

V 

167 

133 

98.89 ± 

3.00 

Ti 

223 

117 

98.45 ± 

2.97 

Sc 

163 

106 

73.64 ± 

2.57 

Ca 

153 

98 

80.32 ± 

2.67 

K 

110 

91 

59.98 ± 

2.31 

Ar 

109 

86 

61.18 ± 

2.32 

Cl 

72 

82 

49.41 ± 

2.09 

s 

80 

78 

59.58 ± 

2.27 

p 

52 

74 

49.82 ± 

2.08 

Si 

50 

71 

72.20 ± 

2.48 

A1 

38 

69 

51.47 ± 

2.10 

Mg 

40 

67 

61.03 ± 

2.27 

Na 

33 

65 

50.17 ± 

2.06 

Ne 

35 

63 

54.55 ± 

2.14 

Lead target 

Mn 

1082 

484 

500.52 ± 

13.42 

Cr 

405 

190 

223.00 ± 

6.18 

V 

250 

160 

130.18 ± 

4.64 

Ti 

335 

142 

135.00 ± 

4.67 

Sc 

244 

129 

104.01 ± 

4.11 

Ca 

230 

120 

98.20 ± 

3.98 

K 

165 

112 

79.76 ± 

3.60 

Ar 

163 

107 

77.23 ± 

3.54 

Cl 

107 

102 

59.97 ± 

3.14 

s 

120 

98 

75.75 ± 

3.47 

p 

78 

94 

63.66 ± 

3.19 

Si 

75 

91 

86.28 ± 

3.65 

A1 

56 

88 

61.90 ± 

3.12 

Mg 

60 

86 

74.14 ± 

3.38 

Na 

49 

84 

66.19 ± 

3.20 


Data from reference 9. 




Table 3. Distribution of Element Production Cross-Section Differences Between 
Theory and Experiment for Beams Fragmenting in Various Targets 

Experimental data used in comparisons are from ref. 11 for 1.88A GeV beams 
and from ref. 9 for 1.55 A GeV beams 


Difference, percent 

Cross sections, percent 

Langley 

NRL 

1.88 A GeV iron beams 

Within error bars 

62 

17 

<25 

15 

8 

26-50 

18 

22 

51-100 

5 

38 

>100 

0 

15 

1.55 A GeV iron beams 

Within error bars 

5 

5 

<10 

19 

2 

11-25 

35 

7 

26-50 

37 

31 

51-100 

5 

34 

>100 

0 

21 


Table 4. Distribution of Elemental and Isotopic Cross-Section Differences 
Between Theory and Experiment for 1.65A GeV Argon Beams 
Fragmenting in Carbon and KC1 Targets 

[Experimental data used in comparisons are from ref. 10] 


Difference, percent 

Cross sections, percent 

Elemental 

Isotopic 

Langley 

NRL 

Langley 

NRL 

Within error bars 

8 

4 

35 

34 

<10 

12 

12 



11-25 

34 

42 

5 

7 

26-50 

27 

15 

13 

18 

51-100 

11 

19 

36 

22 

>100 

8 

8 

11 

19 
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